JA 


OSSIFIED 


405  066 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMtnOH  STAIION  RIEXRNDRU.  VIROINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  puJTxsse 
other  tium  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  -(Aiatsoeverj  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  ri^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


NEW  METHODS  FOR  SYNTHESIS  OF 
NONliNIFORiMLY  SPACED  ANTENNA  ARRAYS 

TECHNICAL  DOCUMENTARY  RERORT  NO.  RADC.TDR.63.131 

May  1963 


Electronic  Warfare  Laboratory 
Rome  Air  Development  Center 
Reieorch  ond  Technology  DIvlilon 
Air  force  8y«temi  Commend 
Orlfflii  Air  force  late,  Now  York 


frolect  No.  46^0,  Totk  No.  454001 


{Propor^fl  under  Contraef  Ne.  Af 30(602)-271 1  by  the  Antenno  Loborofory, 
Disporfmenf  of  Electrical  Erifineerlne,  Ohio  Stote  University  Research 
o'lndotlon,  Celumbue,  Ohio.  Author:  H.  Uni) 


When  Government  drawings,  specifications,  or  other  data  are 
used  for  any  purpose  other  than  in  connection  with  a  definitely 
related  Government  procurement  operation,  the  United  States 
Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever,  and  the  fact  that  the  Goveriunent  may  have  formulated, 
furnished,  or  in  any  way  supplied  the  said  drawings,  specifications, 
or  other  data,  is  not  to  be  regarded  by  implication  or  otherwise  as 
in  any  manner  licensing  the  holder  or  any  other  person  or  corpora^- 
tion,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or 
sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

The  Government  has  the  right  to  reproduce,  use,  and  distribute 
this  report  for  governmental  purposes  in  accordance  with  the  contract 
under  which  the  report  was  produced.  To  protect  the  proprietary 
interests  of  the  contractor  and  to  avoid  jeopardy  of  its  obligations 
to  the  Government,  the  report  may  not  be  released  for  non-governmental 
use  such  as  might  constitute  general  publication  without  the  express 
prior  consent  of  The  Ohio  State  University  llesearch  Foundation. 

Qualified  requesters  may  obtain  copies  of  this  report  from  the 
ASTIA  Document  Service  Center,  Arlington  Hall  Station,  Arlington  12,  ^ 

Virginia.  Department  of  Defense  contractors  must  be  established  for 
ASTIA  services,  or  have  their  "need-to-know"  certified  by  the  cogni¬ 
zant  military  agency  of  their  project  or  contract. 


I  " 


Do  not  return  this  copy.  Retain  or  destroy. 


FOREWORD 


The  reseaxch  reported  herein  was  performed  in  the  Antenna 
Laboratory,  Department  of  Electrical  Engineering,  The  Ohio  State 
University,  under  Contract  No.  AF50(602)“2711  between  The  Ohio 
State  University  Research  Foundation  and  Rome  Air  Development  Center, 
Air  Force  Systems  Command,  USAF,  Griff iss  AFB,  N.  Y.  Mr,  Karl  Kirk 
of  Rome  Air  Development  Center  was  the  contract  initiator. 

This  report  (Contractor's  No.  1425-1)  was  written  by  Dr.  H.  Unz 
of  The  Ohio  State  University  Antenna  Laboratory. 


Much  of  the  research  work  reported  in  this  paper  was  done  while 
the  author  was  a  consultant  for  the  Research  and  Development  Depart¬ 
ment,  Apparatus  Division,  Texas  Instruments,  Incorporated,  Dallas, 
Texas,  and  it  is  published  with  their  permission.  The  author  is  es¬ 
pecially  grateful  to  Mr.  John  B.  Travis  from  Texas  Instruments,  Inc,  , 
for  initiating  this  research  program  and  for  his  support  and  help. 


M  &  W  Jnc»,  Syractae,  New  Yorli 
102  -  5714/63 


ABSTRACT 


Four  new  methods  of  synthesis  of  nonuniformly  spaced  antenna 
arrays  ar^  given; 

L  The  mechanical  quadratures  method,  developed  recently  by 
Bruce  and  Unz  and  independently  by  Lp. 

11.  The  eigenvalues  method. 

III.  The  expansion  method. 

IV.  The  orthogonalization  method  using  the  Schmidt's  procedure. 

An  exponential-decay  directive  pattern  is  suggested  in  order  to 
avoid  numerical  integrations.  A  summary  of  the  research  work  in  non¬ 
uniformly  spaced  antenna  arrays  is  given. 
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NEW  METHODS  FOR  SYNTHESIS  OF 
NONUNIFORMLY  SPACED  ANTENNA  ARRAYS 

1,  INTRODUCTION 

The  importance  of  directive  antennas  was  realized  in  the  early  days 
of  radio  communications.  The  principles  of  wave  interference,  on 
which  systems  of  directive  radio  are  based,  has  been  known  probably 
for  several  centuries.  However,  the  first  thorough  treatment  of  this 
subject  was  conducted  by  Huygens  and  by  Fresnel,  who  established  the 
wave  theory  of  light  in  the  early  part  of  the  nineteenth  century. 

During  the  decade  1920-1930,  a  concentrated  effort  on  directive 
properties  of  antenna  arrays  was  started.  During  this  period  short¬ 
wave  'adio  communications  were  becoming  more  popular,  taking  the 
place  of  long  radio  waves,  and  the  use  of  antenna  arrays  of  reasonable 
size  became  more  feasible. 

In  1937  Wolff^  published  his  method  of  synthesizing,  any  arbitrary 
far-zone  circularly  symmetric  pattern  with  radiators  uniformly  distribu>» 
ted  along  an  array  axis.  His  theory  was  based  upon  comparison  of  the 
far-zone  field  of  a  pair  of  radiators  to  a  term  of  the  Fourier  series  ex¬ 
pansion  of  the  prescribed  pattern. 

During  the  second  world  war,  the  invention  and  the  use  of  radar 
increased  the  interest  in  directive  arrays.  In  1943  Schelkunoff^  util¬ 
ized  the  correspondence  between  nulls  of  the  pattern  of  a  linear  array 
having  equidistant  elements  and  the  roots  of  a  polynomial  in  the  complex 
plane.  Different  types  of  pattern  variations  were  derived  by  choosing 
different  distributions  of  the  zeros  of  the  polynomial. 

In  1946  Dolph^  devised  a  method  of  synthesizing  an  optimum  pattern 
of  isotropic  elements  equally  spaced  in  a  uniform  broadside  linear  array, 
in  which  the  elements  are  fed  in  phase  and  are  symmetrically  arranged 
about  the  center  of  the  array.  The  resultant  current  distribution  across 
the  array  is  based  upon  the  properties  of  the  Tchebyscheff  polynomials 
and  offers,  from  the  design  standpoint,  much  greater  control  of  the 
pattern  than  previous  theories. 


1 


Woodward*  and  Lawson*  gave  a  method  for  calculating  the  field 
over  a  plane  aperture  to  produce  a  given  polar  diagram  and  discussed 
the  theoretical  precisioii  with  which  an  arbitrary  radiation  pattern  may 
be  obtained  from  an  array  of  finite  size.  This  method  is  baaed  on  the 
idea  of  matching  the  radiation  pattern  of  the  aperture  at  a  finite  number 
of  directions  with  the  preecribed  radiation  pattern. 

Van  der  Maas*  gave  the  ideal  Tchebyscheff  space  factor  corres¬ 
ponding  to  a  continuous  source  of  finite  length,  in  the  following  form: 


(la)  F:(0)  =  cos  (Tr  u  u^  -  A*  ) 


2a 

(lb)  u  =  sin  0 

X 

where  2a  in  the  length  of  the  line  source,  \  is  the  wavelength,  and  0  is 
the  angle  from  the  normal  to  the  antenna  axis.  From  the  above  it  can 
be  found  that: 

(lc)  ch(irA)  =  side  lobe  ratio. 


By  using  Fourier  integrals,  Taylor'^  showed  that  the  aperture  current 
distribution  required  in  order  to  produce  the  ideal  radiation  pattern  iii 
Eq.  (1)  is: 


(2) 


irA* 

^:(p)  =  — 


+  2  6(p  -  tt  )  +-2  6(p  +  tt)  ,  p*  Tr‘ 


g(p)  =  0  ,  p*>Tr* 


2 


77  X' 

where  p  =  ■— —  ,,  5  is  the-Dijcaic  delta  function,  and  I  is  the  modified 
a  1 

Bessel  function.  This  pattern  was  considered  later  by  Taylor®  from 

the  complex  function  theory  point  of  view.  Because  of  theoretical  Hmi- 

tations,  the  above  pattern  cannot  be  obtained  from  a  physically  realissable 

antenna^but  its  ideal  characteristics  can  be  approached  arbitrarily 

closely. 

Cheng  and  Ma®^  proposed  a  new  approach  to  the  uniformly  spaced 
linear-array  analysis.  They  considered  the  current  distribution  in 
the  discrete  elements  of  the  linear  array  as  the  sampled  data  values  of 
a  continuous  function,  and  used  known  relations  in  Z  transforms  devel¬ 
oped  for  sampled-data  systeins  in  order  to  express  the  array  poly¬ 
nomial  in  a  closed  form,  where  the  array  properties  are  easily  found. 

In  all  of  the  above  cases,  the:  spacings  of  the  array  elements  are 
assumed  to  be  uniform.  In  1956  Unz^’ ^  ^  ^  first  introduced  linear 

arrays  with  arbitrarily  distributed  elements  and  developed  a  theory  for 
them,  A  relationship  between  the  currents  in  the  elements  of  the  array, 
their  distribution  along  the  axis:  of  the  array,  and  the  coefficients  of  the 
complex  Fourier  expansion  of  the  radiation  pattern  were  given  in  a 
matrix  form  which  involved  Bessel  functions.  It  was  pointed  out  that  a 
nonuniformly  spaced  antenna  array  has  more  degrees  of  freedom  than 
a  similar  uniform  array  with  equally  spaced  elements,  and  as  a  result 
its  performance  should  be  better. 

Further  numerical  work  on  the  subject  has  been  done  by  King, 
Packard  and  Thomas*  ^  ,  who  showed  that  nonuniform  antenna  arrays 
with  unequally  spaced  elements  are  also  more  broadband  for  different 
source  frequencies.  This  is  espeqially  useful  because  of  the  earlier 
development  of  frequency-independent  antenna  elements  by  Rumsey, 
DuHamel,  Isbell,  and  others*®.  Sandler**  showed  some  equivalences 
between  e'  ’ally  and  unequally  spaced  arrays. 

Andreasen*  ®  and  several  of  his  associates  did  extensive  numerical 
work  on  digital  and  analog  computers  in  order  to  find  the  general  be¬ 
havior  of  the  radiation  patterns  of  nonuniformly  spaced  arrays,  in  par¬ 
ticular  when  the  average  spacing  is  larger  than  one  wa  /elength.  Swenson 
and  Lo*  ®  considered  the  use  of  nonuniformly  spaced  arrays  for  large 
radio  telescopes.  Harrington*^  used  perturbational  procedures  for 
reducing  the  side-lobe  level  of  a  nonuniformly  spaced  array  with  uni¬ 
form  excitation.  Extensions  of  the  previous  theory  to  dipole  elements 
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in  the  Fresnel  zone^  ®,  and  to  nonuniform  spacings  larger  than  one 
wavelength^  were  given  by  Unz,  Bruce  and  Unz®  °  gave  a  possible 
condition  for  broadbanding.  Maffett®  ^  discussed  array  factors  with 
a  nonuniform  spacing  parameter  by  using  well-known  numerical  inte¬ 
gration  techniques  (trapezoidal,  Simpson},  and  considered  a  statistical 
theory.  Yen  and  Chow®  ^  ^  discussed  the  possibility  of  expressing  the 
radiation  pattern  of  large  nonuniformly  spaced  arrays  in  closed  form 
by  using  integration  by  stationary  phase  technique, 

Pokrovskii®  ®  designed  a  4r-  and  6-  element  array  having  nonuni¬ 
form  spacing  and  uniform  excitation  of  the  elements,  with  improved 
patterns  over  the  Dolph-Tchebyscheff  array  for  the  same  length  and  the 
same  number  of  elements..  While  his  method  of  solution  is  general, 
the  transcendental  equations  become  more  and  more  involved  for  larger 
arrays.  Additional  discussion  has  been  given  by  Brown®*. 

Nonuniformly  spaced  arrays  have  several  advantages  over  uni¬ 
formly  spaced  arrays;  namely,  their  performance  can  be  better,  fewer 
elements  can  be  used,  and  they  are  more  broadband.  However,  non¬ 
uniformly  spaced  arrays  are  seldom  used  in  practice  at  present.  Be¬ 
sides  being  latecomers,  the  great  difficulty  in  synthesizing  sucl)  arrays 
seem  to  be  the  main  reason  for  the  hesitancy  to,  use  them. 

Recently  a  new  method  of  synthesizing  radiation  patterns  using 
nonuniformly  spaced  arrays  has  been  developed  by  Bruce  and  Unz®  ^ 
They  applied  a  mechanical  quadratures  formulation  in  order  to  trans¬ 
form  a  continuous  aperture  distribution  to  a  nonuniformly  spaced  array. 
Because  of  its  importance,  a  summary  of  this  work  will  be  given  in 
Section!  3  of  this  paper.  Some  work  along  similar  lines  has  been  done 
independently  by  Lo®  *  . 

The  aim  of  this  paper  is  to  give  new  methods  for  synthesis  of  non¬ 
uniformly  spaced  antenna  arrays.  It  is  hoped  that  these  new  methods 
will  be  found  useful  in  the  design  of  nonuniformly  spaced  arrays  and 
helpful  in  achieving  the  general  use  of  their  full  potentialities , 
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2.  THE  EXPONENTIAL  PATTERN 

One  of  the  possible  directive  radiation  patterns  with  exponential 
decay,  which  will  be  used  later  on,  can  be  written  in  the  form: 

-a*  sin*  8 

(3)  .  F(0)  =  e  cos  (2b  sin  0) 

where  0  is  the  angle  of  the  radius  vector  from  the  normal  to  the  array 
axis,  and  "a",  "b"  are  arbitrary  constants.  The  maximum  value  of  the 
pattern  in  the  direction  normal  to  the  array  axis  is  normalized,  F(0)  =  1 . 

The  fir  st  null  of  the  radiation  pattern  will  be  at  0  =  0^  such  that 

IT 

(4a)  2b  sin  ^ 


and  the  beamwidth  w  between  the  first  nulls  Will  be  w^  =  2  6^:  one 

o  o  o’ 

obtains,  therefore. 


(4b) 


4  sin  2 


Using  Eq.  (4b),  the  constant  "b"  can  be  determined  if  the  beamwidth 
between  the  nulls  w^  is  specified. 

Assuming  the  first  side-^lobe  maximum  of  the  radiation  pattern  to 
be  at  0  =  01 ,  then  the  side^lobe  level  M  of  the  radiation  pattern  will  be 
defined  by: 


(5a) 


M 


F(0) 

'  F(0l)  ' 


,  a*  sin*  0, 

1  _  e  * 

F(0i )  cos  (2b  sin  01 ) 


The  side-lobe  level  in  decibels 
(5b)  Mjjjj  =  20  log  M  = 

=  20  [a*  sin*  0j  log  e  -  log 


will  bi;.  then: 


cos  (2b  sin  0i)  ]  . 


5 


In  order  to  find  tlie  position,  of  the  fir  st  maximuna,  phe  should  solve 
F'(9i)  =  0.  Using  Eq.  (3)  one  obtains; 

b  sin  (Zb  sin  8i)  +  sinSi  cos  (2b  sin  Gj)  =  0. 


Using  the  notation  a  =  2b  sin  ©i  one  gets: 


(6) 


tg,a  „  ^  a^ 
a  2b? 


Using  Eq.  (6)  and  the  same  notation,  Eq. 

'^db  I 

— =  “  0.217  a  tg  a  -  log  I  cos  a 


(5b)  becomes: 


Equations  (6)  and  (7)  may  be  used  in  order  to  find  the  relationship 
between  the  constants  "a",  "b",  and  the  side-lobe  level  and  the 

results  are  shown  graphically  in  Figures  1,  2,  3. 

The  required  radiation  pattern  is  given  usually  in  terms  of  its 
beamwidth  w  and  the  side-lobe  level  From  Eq.  (4b)  the  constant 

"b"  of  the  raHiation  pattern  can  be  found  from  the  beamwidth  w^.  For 
a  given  side-lobe  level  M^|^,the  transcendental  equation  (7)  can  be 
solved  for  a.  This  could  be  accomplished  simply  by  using  Figure  2. 

Then  Eq.  (6)  and  the  value  of  the  constant  "b”  should  be  used  in  order 
to  find  the  constant  "a”.  Figure  1  could  be  helpful  in  this  case.  Figure 
3,  which  is  a  combination  of  the  previous  figures,  gives  a  direct  rela¬ 
tionship,  in  a  graphical  form,,  between  the  side-lobe  level  and  the 

a^ 

constant  . 

We  have  seen  that  if  the  beamwidth  w^  and  the  side-lobe  level 
are  given,  the  constants  "a",  "b"  may  be  found  and  the  directive,  ex¬ 
ponential-decay  radiation  pattern  in  Eq.  (3)  will  be  completely  specified. 

3.  MECHANICAL  QUADRATURES  METHOD 

Recently  a  new  method  of  synthesis  of  nonuniformly  spaced  linear 
arrays  has  been  developed  by  Bruce  and  Unz?  *  and  independently 
by  Lo^  *  .  They  applied  mechanical  quadratures*  ^  in  order  to  transform 
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a  continuous  aperture  distribution  to  a  nonUniformly  spaced  array. 
Beicause  of  its  importance,  a  sumniary  of  this  new  method  ^  3  ,  2  4  j_g 
presented  here. 

Using  weight  functions ,,  the  general  form  of  tneehanical  quadra¬ 
tures*  ^  is  given  by: 

+.1.  n 

(8) 


w(x)  f  (x)  dx 


The  n  values  of  aj  are  real,  distinct,  within  the  interval  (-1,  +1),  and 

are  the  roots  of  a  polynomial  orthogonal  with  respect  to  the  weight 
function  w(x)  on  the  interval  (-1,  +1).  The  n  values  of  Hj  are  real, 
positive,  and  can  be  determined  by  solving  the  set  of  linear  equations 
generated  by  letting  f(x)  =  x^  for  k  =  0,  1,  2  •••  2n  -  1.  The  only  re¬ 
quirements  on  the  weight  function  are  that  it  does  not  vanish  within  the 
interval  of  integration  and  that  it  is  integrable. 

For  a  finite  aperture  with  symmetric  excitation,  the  radiation 
pattern  can  be  written  as: 

+1 

(9)  F(u)  =  ^  g(x)  cos  uxdx 

-1 


where  the  aperture  length  has  been  normalized,  u  =  it  sin 6,  and  x  is 
the  distance  from  the  origin  measured  in  half  wavelengths. 

Comparing  Eqs.  (8)  and  (9),  and  taking  the  weight  function  to  be 
equal  to  the  aperture  distribution  w(x)  =  g(x),  one  obtains: 

n 

(10)  F(u)  =  ^  HjC08(uap. 

j=l 


Equation  (10)  is  identified  as  the  nonuniform  symmetric  array  radiation 
pattern,  where  the  element  positions  are  given  by  the  n  values  of  a.  and 

the  current  excitation  coefficients  are  given  by  n  values  of  H., 

J 
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1 

'\ 

1 


(11) 


In  case  of  no  weight  fuhctibh,  w(x)  =  1,  Eq,  (8)  becomes: 
+1  h 

^  f(x)dx  =  ^  Hjif(aj):  ‘ 

j=l 


-1 


and  as  a  result,  Eq.  (10)  becomes: 


n 

(12)  F(u)  =  ^  Hj  g(aj)  cos  (ua-j)  . 

j=l 


Bruce  and  Unz^  3,2  4  applied  this  method  to  the  exponential  direc¬ 
tive  radiation  pattern  described  in  Section  2,  as  well  as  to  the  Dolph- 
Tchebyscheff  pattern^with  very  good  results.  In  the  case  of  the^ ex¬ 
ponential  pattern,  by  taking  the  weight  function  to  be  w(x)  =  e“^  one 

finds  the  values  of  a.  and  Hi  to  be  tabulated^  ^  . 

J  J 


The  basic  requirement  for  the  mechanical  quadratures  method  is 
the  knowledge  of  the  aperture  distribution  g(x)  in  Eq.  (9).  If  the  radia¬ 
tion  pattern  F(u)  is  prescribed,  g(x)  could  be  found  sometimes  by  using 
the  inverse  Fourier  transform,  as  for  the  case  of  the  exponential  pat¬ 
tern^  3,2  4  Section  2.  Sometimes  the  aperture  distribution  g(x)  may 
be  found  by  taking  the  envelope^ of  an  already  designed  uniformly- 
spaced  array,  as  for  the  case  of  the  Dolph-Tchebyscheff  pattern^ 

In  both  case's  the  current  at  the  aperture  edge  will  give  difficulty,  as 
can  be  seen  in  Eq.  (2),  for  example,  and  should  be  discussed  separ¬ 
ately*  ,  “*  * . 
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For  the  design  of  a  general  prescribed  radiation  pattern,  the  diffi¬ 
culties  of  this  method  are  twofold;  (a)  To  find  the  aperture  excitation 
from  the  prescribed  radiation  pattern  by  a  Fourier  integral;  (b)  Approsd- 
mating  the  aperture  excitation  by  a  polynomial,  a  process  which  will 
involve  the  inversion  of  a  matrix.  However,  for  directive  patterns 
certain  short  cuts  in  the  process  have  been  shown  above ^ 
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4.  EIGENVALUES  METHOD 

A  symmetric ,  nbnuniformly  spaced  antbnna  array  will  give  the 
radiation  pattern: 


(13) 


F(u) 


^  A|  cos  (uX|) 


l~0 


where  u  =  j^sinQ,  6  being  the  angle  with  the;  normal  to  the  array  axis, 
and  x^  is  the  disttince  from  the  center,,  measured  in  half  wavelengths. 
In  Eq.  (13)  the  radiation  pattern  F(u)  is  given  and  the  array  elenient 
distribution  X|  and  the  element  amplitudes.  A|  have  to  be  deterrnined. 

Let  us  discuss  the  integral: 

-hr 

{14a)  I(X|;x^) 

-7r 


=  y  cos  (ux|)  cos  (ux^)  du. 


When  X| 


=  Xjjj  one  obtains: 


sin  2  x^  7T 

(14b)  I^  =  I(x^;x^)  =  7r[  1  +  2x^7r  ^  ’ 


By  using  trigonometric  identities  Eq.  {14a)  becomes: 
sin(x|  +Xjj^)7r  sin(x^ 


X  Itt 
m' 


I(X|:Xm)  = 


^ 


Equating  the  last  expression  to  zero  leads  to: 
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Rearranging, 

Xj{  [  sin  (x^  +  x^)  7T  +  sin  Cx^  _-  x^)  tt]  = 

=  x^  [sin  (x^^  +  x^)7r  -  sin  (x^'  -  x^)  tt]  . 


Using  trigonometric  identities,  one  obtains  finally: 


(1  4c)  (x^  7r)  tg  (x^  7r)  =  tg  (x^ir)  . 


It  is  found  that  for  the  eigenvalues  of  Eq.  (14c),  the  integral  I  (x,  ;  x  ) 

I  m' 

in  Eq.  (14a)  will  have  orthogonality  properties  such  that; 


(14d)  I(Xjj  ;  x^)  =  7t[  1  + 


sin  2  X , 


7r 


-1  6 


i  ;m 


where 


1  i  =  m 

^0  f  Km 


is  the  Kronecker  delta. 


From  Eq.  ()''»,<  for  the  eigenvalues  one  can  see  that  all  the  posi¬ 
tions  of  the  radiating  elements  are  determined  by  the  position  of  the 
first  element  with  respect  to  the  middle  of  the  array.  In  order  to  in¬ 
clude  the  first  eigenfunction  in  the  complete  set,  the  first  element 
position  should  be  0  _<  x^  <  1  ;  in  other  words,  the  distance  between 
the  two  innermost  elements  should  be  smaller  than  one  wavelength.  The 
eigenvalues  in  general  are  unequally  spaced  for  small  values  and  are 
almost  equally  spaced,  a  half  wavelength  apart,  for  large  values. 

Using  the  orthogonality, properties  in  Eqs.  (14)  one  may  obtain 
from  Eq.  (1  3): 


(15) 


I  TT  *• 


sin  2  X.  TT  -1  p  +7T 

- ]  \  F(u)  cos  (ux  )  du  . 

2x  7r  J  ' 

*  “7T 
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Thus,  one  could  get  a  nonunifbrmly  spaced  array  which  will  produce 
any  required  radiation  pattern..  Its  element  distribution  will  be  deter¬ 
mined  by  Eq.  (14c)  and  the  currents  will  be  determined  by  Eq.  (15). 

For  the  particular  case  of  the  exponential  pattern,  Eq.  (3)  may  be 
rewritten: 


-T  u.  ,, 

(16)  F(u)=e  cos  ( —  u)  f 

Tt' 

1,  'i' 

where  u  =  tt  sin©.  Substituting  Eq.  (16)  into  Eq.  (15)  and  using  trigo¬ 
nometric  identities,  one  may  obtain  definite  integrals,  which  may  be 
evaluated  approximately^’  for  a  >  2: 


TT  2 

P  -  —2“ 

^1^)  \  6)  Cos  (2pu) 

*^0 


^  ^  e  ^  cos(2pu)du=  e  ^ 

0 


Using  Eq.  (17)  one  may  obtain  for  the  exponential  pattern: 


_  sin  2  X  71 

A,  =  [  1  +  1  -1  [ 

I  2  a  f  ^  ■*  ^ 


b  +2  X 


(18a) 


i\Z 


K  2 
b  -  2  X 


r  -  (- 

+  e. 


]  . 


or  in  alternative  form: 


(1  8b) 


o 

2  a 


[1  + 


2x.  7r 


n...  2  x^77  /  ■^vz  , 

J  e  e  ch  — ^)  . 
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Equations  (18)  give  the  current  distribution  of  an  antenna  array  which 
produces  the  exponential  radiation  pattern  ih  Eq.  (3)  br  (16),  and  the 
elements  of  the  array  are  distfibuted  according  tP  the  eigenvalues  of 
Eq.  (I4c), 

From  Eq.,  (18a)  one  can  see  that  for 
( 19)  Y  ^  31  +  b 

the  value  of  quickly  becomes  very  small,  and  the  contributions  of 
these  terms  to  the  radiation  pattern  are  so  small  that  they  could  be 
neglected.  Equation  (19)  gives  an  . indication  of  the  number  of  elements 
required  to  produce  the  exponential  pattern  within  a  certain  approxi¬ 
mation. 

For  the  particular  condition  of  uniformly  spaced  arrays ,  one  could 
distinguish  between  two  cases: 

A.  Odd  number  of  elements: 


In  this  case  x^  =  0  and  from  Eq.  (I4c)  one  obtains  x^  =  f .  Equation 
(  15)  becomes  then: 


(20a) 


.  +  7T 


L  =.  —  r  F(u)  cos  (i  u)  du 
i  Tt  J 

"  TT 


t  >  0 


(20b) 


1 

V  =  F(»)du  .. 

-  TT 


B..  Even  number  of  elements; 

In  this  case  x^  = and  from  Eq.  (14c)  one  obtains  x^  =  f  + -^ 
Equation  (15)  becomes  then: 


(21) 


1  c'*"^  1 

I  =  ~  \  F(u)cos(i  +-^)  udu 


TT 
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«“atA«vC'i«sKs^ 


and  the  relationships  for  the  exponential  pattern  should  be  changed 
accordingly,  in  Eqs.  (18)., 

In  the  last  two  methods  of  Sections  3  and  4  the  positions  of  the 
array  elements  have  been  determined  during  the  synthesis  procedures. 
These  methods  will  fail  in  general  when  the  element  positions  are  pre¬ 
assigned.  In  the  next  two  sections  we  will  discuss  cases  of  this  kind. 

5.  EXPANSION  METHOD 

When  the  radiation  pattern  F  {xn)  is  given  and  the  nonuniformly  dis¬ 
tributed  antenna  array  element  positions  are  preassigned,  the  determi¬ 
nation  of  the  currents  in  Eq.  (13)  is  relatively  difficult.  The  reason 
is  that  one  is  required  to  expand  a  given  function  in  terms,  of  a  non- 
prthogonal  set  of  functions.  This  problem  has  been  discussed  by  Kan-r 
tbrovich  and  Krylov^ 

One  possible  method  of  evaluating  the  coefficients  A^  in  Eq.  (13) 
will  be  to  take  any  complete  set  of  functions  {V'^('^)}>  naultiply  both 

sides  of  Eq.  (13)  by  ^j^(u)  and  integrate.  One  will  obtain  then  a  system 
of  equations; 

L 

(22a)  /  cn  -  A«  =  b  '  m  -  0,  17  2‘''’L 

Z..  f;m  *  m  ’  ’ 

j  =  0 

where: 

r+'/r 

(22b)  =  \  cos  (ux^)  (u)  du 

~  7r 


( 

(22c)  b^=  F(n)f^{u)du 

“  7T 

where  c.  and  b  are  known  and  A,  is  to  be  evaluated  by  solving 
£^m  m  .  I  /  s 

(L  +  1)  linear  equations  with  L  +  1  unknowns. 
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It  will  be  useful  to  indicate  also  another  method  of  obtaining  sys¬ 
tem  (Z2a),  which  is  doubtless  even  more  convenient.  Let  us  expand 
Eq.  (13)  in  terms  of  any  orthonormal  system  Let: 

00 

(23a)  F(u)=J 

m=0 


(23b)  cos(ux^)=^ 

m=0 


where  b  is  given  by  Eq.  (22c)  and  c.  is  given  by  Eq.  (22b).  Sub- 

m  X  j  in 

stituting  Eqs.  (23)  into  Eq.  (13)  and  equating  coefficients  of'^u)  in  both 
sides,  one  obtains  Eq.  (22a).  We  note  that  the  second  method  of  ob¬ 
taining  system  (22a)  presupposes,  as  distinct  from  the  first,  a  system 
of  functions  i^{^)  that  is  orthogonal  and  normalized. 

One  possible  such  set  is  -  cos  mu.  This  set  will  have  a 

definite  advantage,  as  Eq.  (22a)  will  be  simplified  considerably  if  some 
of  the  elements  of  the  array  are  uniformly  distributed  at  multiples  of  a 
half  wavelength, apart.  Also,  the  integrals  (22b)  are  readily  available. 

The  above  method  could  be  extended  to  the  case  of  asymmetric 
nonuniformly  spaced  arrays: 

V 

(24)  F(u)  =  2  e  • 

^  =  0 

One  could  take  then  ^j^(u)  =  Repeating  the  above  process,  one 

will  obtain: 


(25a) 


no  =  0,  1,  2  L 


(25b) 


m 


imu  sin7r(x. +m) 

e  du :=  2  ir - 

7T  +m) 


i 

1 
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,  n  +Tr  .  „ 

(25e)b^=.|  F(u)  e^'^du 

'  "TT 


In  this  case,  in  general,  b^and  A  .  will  be  complex  and  one  should  take 
both  +  (m)  and  -  (m).  The  theory  developed  by  Unz^’*  ^  ^  originally  is  a 
particular  case  of  the  method  described  in  this  section, 

6,  ORTHOGONALIZATION  METHOD 

in  order  to  expand  an  arbitrary  function  in  terms  of  a  nonorthogonai  . 
set  of  functions,  the  orthogonalization  procedure  by  Schrnidt  may  be 
used.  This  procedure  has  been  discussed  in  detail  by  Kantorovich  and 
Krylov^®,  and  its  details  for  complex  functions  may  be  found  in  Appen¬ 
dix  A  of  the  present  paper. 

It  may  be  shown^°  that  from  an  infinite  set  of  nonorthogonai  complex 
functions: 

<t’i  (n),  $2  (u)  ••• 

an  orthonor'mal  set  of  functions 


t 


! 


Tpi  (u),  fz  (n)  ••• 

may  be  derived  by  using  the  relationships: 


(26a)  fi  (u)  = 


<t'i  (u) 


\  <)),  (u)  (t)f(u)  du 

r  a 


{26b) 


n  ^ 

‘^n+l  ^n+1 

j  =  l-  J 


ff’’! 

■  ,  1 

•j  ' 

a  / 

^  L — 

j  =  l 


where  AA*  =  |  a|^  .  The  orthonormal  functions  if/^  (u), 
may  be  found  in  succession. 
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The  brthonormal  set  ^^(u)  obeys: 


(27>  j" 


m 


1  I  =  m 


.0  i  ^ 


m 


Let  UiS  assume  that  a  set  of  nonorthogdnai  functions  is  given  by: 


(28)  4>^tv.)  e^^nV 


n  =  1,  2,  •"  L 


where  Xj  <  <  X3  <  ••*  <  x^.  We  would  like  to  orthogonalize  this 

set  of  functions,  in  the  region  {-71-,  +77),  where  the  orthonor.maljs'nt.Will 
be  denoted  by  V'nCu)-  From  (28)  and  (26a)  one  obtains  : 


{29a)  V'i(’^) 


i _  " 


nJ2 


TT 


Let  us  define  the  following  function: 

(29b)  f  ^ 

‘^-TT 


+  7T 


"  7T 


8in{x^  -  X^)7T 

-  2.77  - -  -  27rS|na  - 

{Xi  -  X^)  7r 

where  £  -  1. 

Using  Eqs.  (28),  (29),  and  Eq,*  (26b), foT' n  =  I  ,  {ohe} obtains:  . 


(  29c)  '  fz,  (u)  = 


-  Si  2  e 
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Similarly,  one  may  find  in  Eq.  (26b)  after  reaffanging: 


(  29<^  ^3  (u)  = 

1  (1  -  Si^2  )e^  -  (S,  z  -  St,  2  Si  3  )e"  -  feSi  Si  aS  ^  3  )  e^  ^ 

^  [I  ^Sf2  ]  ^  [i  -S  L  -S  h  -4  3  +  2Si  3  p 

By  continuing:  the  above  process  in  Eq.  {26b)  one  may  find  forfriaily ; 


(29e) 

f4  (U) 

u 

II 

*  ’‘I "  +  c,  '''•  «*  *2  “  +  c,  >e*  *>  "  +  c.***a‘ 

or  in 

general: 

n 

(29f) 

=  I 

(n)  ix<u 
f. 

<=i 

where  the  constants  c^^”^  may  be  found  by  the  Schmidt  orthogonalization 
procedure  describedabove. 


In  the  case  of  a  general,  asymmetric  nonuniformly  spaced  array, 
one  has; 

L 


(  30  )  F(u) 


=.  I 

£=1 


Taking  the  orthogonal  set  of  functions  ^^^(u)  described  in  Eqs.  (  29  )  one 
could  expand: 

N 

(31a)  F(u)  =  Z  B„V^^(u) 
n=l 
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where,  using  the  orthogonality  relationship  in  Eq.  (27),  one  can  find: 
^  +7r 

(31b)  ®n~  )  F(u)  ijt'*  (u)  d  u  . 

Since  F(u)  is  prescribed,  B  may  be  evaluated  from  Eq.  (31b)  using 
Eqs.  (29).  ^ 

Substitviting  Eq.  (29f),  into  Eq.  (31a)  one  obtains,  after  rearranging:: 
N  n 

(  32  )  F(u)  =  5  I  • 

n=l  f=l 


Comparing  Eqs,  (  32  )  and  (  30  )  one  obtains: 


(33a) 

Ai  =  Bi  ci‘  ^  +  B2 

^+B3  c!^U*- 

%ci 

(33b) 

Aj  =  B2 

C2 

...  B  c^^) 

(33c) 

A}  = 

+ 

...  B 

(33  d) 

II 

C 

B  c<^> 
N  N 

By  using  the  above  procedure  one  can  see  that  in  principle  the  non- 
uniformly  spaced  array  >  an  be  designed  as  well  as  uniformly  spaced 
arrays  without  using  the  process  of  inversion  of  large  matrices.  The 
process  of  finding  the  orthogonal  set  of  functions  in  general  is  rather 
involved  algebraically,  but  they  could  be  found  with  the  help  of  a  digital 
computer  using  the  successive  process  described  by  Eqs.  (26).  In  gen¬ 
eral  the  result  will  correspond  to  the  radiation  pattern  within  a  certain 
approximation,  rather  than  be  equal  to  it. 

The  great  importance  of  this  process  of  solution  of  nonuniformly 
spaced  arrays  is  that  it  is  especially  designated  for  asymmetric  non¬ 
uniformly  spaced  arrays,  which  have  more  degrees  of  freedom  than 
symmetric  nonuniformly  spaced  arrays.  Of  course,  symmetric  arrays 
also  can  be  designed  by  using  this  method. 
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For  the  case  of  an  exponential  pattern  described  in  Section  Z,  one 
is  able  to  use  Eqs.  (l6)  and  (17)  in  order  to  integrate  in  Eq.  (:31b)  and 
find  the  corresponding  coefficients  ^ 

For  the  particular  case  of  uniformly  spaced  arrays  with  y  spacings, 
one  will  have  in  Eq.  (  29b)  : 


(34a) 

X-  -  X  =  f  -  m  ; 

£  m  £m 

and  Eqs. 

(29)  become; 

(34b) 

1  in 

tiin)  =..  e 

2  7f 

(34c) 

,  .  1  i2u 

^2  (u)  =  e 

■J  ^TT 

(344) 

ih  1  \  1  j  3  u 

vj  Ztt 

{34e) 

/  V  1  4nu 

=  - -  « 

\1 2  7r 

and  Eqs,  (  31  ),  (  32  )  and  (33  )  will  reduce  to  the  standard  Fourier 
series  analysis  of  uniformly  spaced  arrays. 

7 .  SUMMARY 

In  this  paper  the  published  research  vmrk  in  nonuniformly  spaced 
antenna  arrays  has  been  summarized  to  date.  Four  new  methods  for 
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synthesis  of  npnunifbrmly  spaced  antenna  arrays  for  a  given  radiation 
pattern  have  been  suggested: 

I.  The  mechanical  quadratures  method,  developed  recently  by 

Bruce  and  Unz^^  ’  ^  *  and  independently  by  Lo^  *  is  summarized, 

II.  The  eigenvalues  method,  where  the  positions  of  the  elements 
are  given  as  eigenvalues  of  a  transcendental  equation. 

III.  The  e3{pansion  method,  which  involves  the  inversion  of  ma¬ 
trices  in  order  to  find  the  currents  of  the  elements. 

IV.  The  orthogonalization  method  by  Schmidt’s  procedure. 

In  methods  I  and  11  the  positions  of  the  elements  are  predetermined 
by  the  weight  function  involved  (method  I),  or  by  the  position  of  the 
innermost  element  (method  11).  In  methods  III  and  IV  the  array  ele- 
ments  are  first  arbitrarily  distributed  (in  accordance  with  the  pre¬ 
requisites,  e,  g..  ,  broadbanding)  and  then  the  current  distributions  in 
the  elements  are  determined.  Methods  I  and  III  will  involve  in  general 
the  inversion  of  a  matrix,  while  in  methods  II  and  IV  this  is  avoided. 

An  exponential- decay  radiation  pattern  is  suggested  and  used  as  an 
example  in  the  above  methods.  Its  main  advantage  in  the  present  syn¬ 
thesis  is  that  the  definite  integrals'  .nvolved  may  be  approximately 
evaluated  explicitly,  and  thus  simplify  the  numerical  work,.  This 
directive  exponential  pattern  could  also  be  used  for  electronic  scanning*^ 
problems  in  the  form: 


(  35  )  F(0  -  y)  =  e  ^^cos  [  2  b  sin  (6  -  y)  ] 


where  y  =  y{t)  is  the  angle  direction  in  which  the  radiation  pattern  will 
have  a  maximum  and  y(t)  is  a  function  of  time.  The  above  methods 
I  could  be  modified  slightly  in  order  to  calculate  the  required  current 
distribution  by  numerical  integration.  The  linear  phase  delay  case  for 
scanning  may  be  found  as  a  particular  case  under  certain  approximations. 

Some  numerical  work  using  the  above  methods  has  been  done,  and 
it  will  be  published  in  the  future  together  with  additional  numerical,  work 
and  comparisons  between  the  four  methods  of  synthesis,  with  regard  to 
the  best  approximations. 
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Sharjp**  and  Willey®^  suggested  methods  of  design  of  linear  and 
planar  arrays  in  order  to  reduce  the  number  of  the  elements  required 
for  a  specific  radiation  pattern.  Ishimaru’*  has  recently  suggested 
the  use  of  the  Poisson’s  sum  formula  for  the  design  of  nonuniformly 
spaced  arrays.  His  method  is  useful  in  treating  npnuhiform  arrays 
with  large  number  of  elements  and  unequally  spaced  arrays  on  curved 
surfaces. 

It  is  hoped  that  by  using  the  synthesis  methods  suggested  by  the 
author  in  the  present  pa^per  and  others,  the  design  of  the  nonuniformly 
spaced  arrays  will  become  simpler  and  their  use  more  generally 
accepted^ 
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APPENDIX 


In  the  following  we  will  es^end  the  orthbgonalization  pTpeess  of 
Schmidt,  as  giveh  by  Kahtorovich  and  Krylov,  to  a  set  of  complex 
functions. 


Let  there  be  given  an  infinite  system  of  complex  functions. 


(A-I)  4ii{u).  i^zM, 

defined  and  continuous  in  the  interval  (a,  b).  We  can  exclude  from  the 
given  system  of  functions  those  that  represent  linear  combinations  of 
the  preceding  ones,  since  by  their  nature  they  do  not  extend  the  sys¬ 
tem.  Let  us  now  carry  through  the.  orthQgonalization  of  the  system.. 
(A-1)  ,  i.  e.  ,  let  us  now  construct  the  orthonormal  system  of  fimctions 


(A- 2)  1}/^  (u),  (u),  ••• 


such  that 


(A- 3) 


(u)  (u)  d  u 


^f.m 


where  5  is  i'he  Kronecker  delta,  and  f*  (u)  is  the  complex  conju- 

f  ,m  na 

gate  function.  Each  function  of  Eq.  (A-2)  is  to  represent  some  linear 

combination  of  the  functions  of  system  (A-1),  i.e,  ,  ij/  (u)  will  have  the 

form,:  ..  .......  ^ 

(A-4)  ^^u  =  a^j^)  <t)  i(u)  +  a^”^  ifz  (u)  +  •••  a|^"^  4,^(u) 


Let  us  carry  out  this  orthogonalizatipn  step  by  step.  The  first  function 
ij/^  (u)  must  have  the  form  c  (t>^  (u). 

Determining  the  constant  c  from  the  condition 


(.»;)  (u)  du  =  1 
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I..'. 


?  ' 
f  >' 


we  fiiid- 


<t>i(u) 


(A- 5)  y/,(u)=  - 

-pb  -| 

F' 

V  ti  (tt)  4>*  (u)  d  u 

y  a.  ^  _ 

Let  us  assume  that  the  first  n  functions  if/i  (u),  (h)>  *'*  ^  (u)  have 

been  determined.  The  function  j^(u)  must  be  a  linear  combination 
of  these  and  the  function  ^(u)  in  the  form: 

(A-6)  +  +  ■"  ^‘^n^n('^)  +1  - 

We  determine  the  constant  c^^  (i  =  1,  2,  •••  n)  from  the  conditions  of 

orthogonality  of  ^i(u),  (u)  •••  multiplying  Eq.  {A-6) 

by  tj/*  (u)  and  integrating  we  obta-in,  using  the  orthogonality  and  normality 
con 


(A 


"la.)  Cl  +  c  J  du  =  0  . 


Sirnilarly,  by  multiplying  by  f*  (u)  and  integrating: 
(A- 7b)  cj  +  c  J  (u)du  =  0  , 


and  so  on. :  Finally  we  obtain 

(Ar-7c)  +  c  J  '*'n  +  l^'^^  ^n  =  0  . 

Substituting  c^  ,  Cj  ,  •*'  from  Ec^.  (A- 7)  into  Eq.,  (A-6)  one  obtains: 

n 

(A-8)  ^n  +  i(")  =  c  I'n  +  iH  "  f  *n+l  t V'j  ^  • 
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*>• 
s: 
s  ^ 

'I; 

r 


The  constant  c. can  be  deternjiined  from  the  normalization  ^  V'n+1  V'n'-fi  ^ 
Using  this,  the  final  result  will  be: 


n 


{A-9y  = 


n  ,  b 

^n+1"*'”  '^n+l^!'^'^)^^'^) 


du) 


j-t 


where  A  ^  =  AA#. 


By  means  of  Eq.  (A- 9)  the  functions  ^2  (u),  ■^3  (0)  •••  may  be  found 

ill  succession  by  cascade  procedure. 

After  having  constructed  the  orthonormal  system,,  one  can  write, 
for  any  arbitrary  function  f{u),  its  series: 


(A- 10a)  f(u)  ^  :Z  Aj^^^(u) 

n=l 


where  the  coefficients  A^  may  be  found  by  using  Eq.  {A«3)  to  be 


(A-lOb)  ^n  =  J  f(u)V£*(u)du 


It  is,  of  course,  generally  impossible  to  guarantee  the  convergence  of  the 
series  in  (A- 10a).  It  would  be,  therefore,  more  correct  to  say  that  this 
series  corresponds  to  the  function  f(H)  rather  than  is  equal  to  it. 

Substituting  Eq.  (A-4)  into  Eq.  {A- 10a)  one  obtains: 


(A- 11)  f(u)  i:^ 


z_ 

n=l 


a["^  <t>i  (u)  +ai"^ 


f(u)  is  written  in  terms  of  combination  of  the  original  nonorthogonal  set. 
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